Computed tomography-based 3D convolutional neural network deep learning model for predicting micropapillary or solid growth pattern of invasive lung adenocarcinoma.

PURPOSE: To investigate the value of a computed tomography (CT)-based deep learning (DL) model to predict the presence of micropapillary or solid (M/S) growth pattern in invasive lung adenocarcinoma (ILADC). MATERIALS AND METHODS: From June 2019 to October 2022, 617 patients with ILADC who underwent preoperative chest CT scans in our institution were randomly placed into training and internal validation sets in a 4:1 ratio, and 353 patients with ILADC from another institution were included as an external validation set. Then, a self-paced learning (SPL) 3D Net was used to establish two DL models: model 1 was used to predict the M/S growth pattern in ILADC, and model 2 was used to predict that pattern in ≤ 2-cm-diameter ILADC. RESULTS: For model 1, the training cohort's area under the curve (AUC), accuracy, recall, precision, and F1-score were 0.924, 0.845, 0.851, 0.842, and 0.843; the internal validation cohort's were 0.807, 0.744, 0.756, 0.750, and 0.743; and the external validation cohort's were 0.857, 0.805, 0.804, 0.806, and 0.804, respectively. For model 2, the training cohort's AUC, accuracy, recall, precision, and F1-score were 0.946, 0.858, 0.881,0.844, and 0.851; the internal validation cohort's were 0.869, 0.809, 0.786, 0.794, and 0.790; and the external validation cohort's were 0.831, 0.792, 0.789, 0.790, and 0.790, respectively. The SPL 3D Net model performed better than the ResNet34, ResNet50, ResNeXt50, and DenseNet121 models. CONCLUSION: The CT-based DL model performed well as a noninvasive screening tool capable of reliably detecting and distinguishing the subtypes of ILADC, even in small-sized tumors.

Prediction of microvascular invasion and pathological differentiation of hepatocellular carcinoma based on a deep learning model.

PURPOSE: To develop a deep learning (DL) model based on preoperative contrast-enhanced computed tomography (CECT) images to predict microvascular invasion (MVI) and pathological differentiation of hepatocellular carcinoma (HCC). METHODS: This retrospective study included 640 consecutive patients who underwent surgical resection and were pathologically diagnosed with HCC at two medical institutions from April 2017 to May 2022. CECT images and relevant clinical parameters were collected. All the data were divided into 368 training sets, 138 test sets and 134 validation sets. Through DL, a segmentation model was used to obtain a region of interest (ROI) of the liver, and a classification model was established to predict the pathological status of HCC. RESULTS: The liver segmentation model based on the 3D U-Network had a mean intersection over union (mIoU) score of 0.9120 and a Dice score of 0.9473. Among all the classification prediction models based on the Swin transformer, the fusion models combining image information and clinical parameters exhibited the best performance. The area under the curve (AUC) of the fusion model for predicting the MVI status was 0.941, its accuracy was 0.917, and its specificity was 0.908. The AUC values of the fusion model for predicting poorly differentiated, moderately differentiated and highly differentiated HCC based on the test set were 0.962, 0.957 and 0.996, respectively. CONCLUSION: The established DL models established can be used to noninvasively and effectively predict the MVI status and the degree of pathological differentiation of HCC, and aid in clinical diagnosis and treatment.

Deep learning-assisted LI-RADS grading and distinguishing hepatocellular carcinoma (HCC) from non-HCC based on multiphase CT: a two-center study.

OBJECTIVES: To develop a deep learning (DL) method that can determine the Liver Imaging Reporting and Data System (LI-RADS) grading of high-risk liver lesions and distinguish hepatocellular carcinoma (HCC) from non-HCC based on multiphase CT. METHODS: This retrospective study included 1049 patients with 1082 lesions from two independent hospitals that were pathologically confirmed as HCC or non-HCC. All patients underwent a four-phase CT imaging protocol. All lesions were graded (LR 4/5/M) by radiologists and divided into an internal (n = 886) and external cohort (n = 196) based on the examination date. In the internal cohort, Swin-Transformer based on different CT protocols were trained and tested for their ability to LI-RADS grading and distinguish HCC from non-HCC, and then validated in the external cohort. We further developed a combined model with the optimal protocol and clinical information for distinguishing HCC from non-HCC. RESULTS: In the test and external validation cohorts, the three-phase protocol without pre-contrast showed κ values of 0.6094 and 0.4845 for LI-RADS grading, and its accuracy was 0.8371 and 0.8061, while the accuracy of the radiologist was 0.8596 and 0.8622, respectively. The AUCs in distinguishing HCC from non-HCC were 0.865 and 0.715 in the test and external validation cohorts, while those of the combined model were 0.887 and 0.808. CONCLUSION: The Swin-Transformer based on three-phase CT protocol without pre-contrast could feasibly simplify LI-RADS grading and distinguish HCC from non-HCC. Furthermore, the DL model have the potential in accurately distinguishing HCC from non-HCC using imaging and highly characteristic clinical data as inputs. CLINICAL RELEVANCE STATEMENT: The application of deep learning model for multiphase CT has proven to improve the clinical applicability of the Liver Imaging Reporting and Data System and provide support to optimize the management of patients with liver diseases. KEY POINTS: • Deep learning (DL) simplifies LI-RADS grading and helps distinguish hepatocellular carcinoma (HCC) from non-HCC. • The Swin-Transformer based on the three-phase CT protocol without pre-contrast outperformed other CT protocols. • The Swin-Transformer provide help in distinguishing HCC from non-HCC by using CT and characteristic clinical information as inputs.

Non-contrast CT radiomics and machine learning for outcomes prediction of patients with acute ischemic stroke receiving conventional treatment.

PURPOSE: Accurate prediction of outcomes for patients with acute ischemic stroke (AIS) is crucial for clinical decision-making. In this study, we developed prediction models based on non-contrast computed tomography (NCCT) radiomics and clinical features to predict the modified Rankin Scale (mRS) six months after hospital discharge. METHOD: A two-center retrospective cohort of 240 AIS patients receiving conventional treatment was included. Radiomics features of the infarct area were extracted from baseline NCCT scans. We applied Kruskal-Wallis (KW) test and recursive feature elimination (RFE) to select features for developing clinical, radiomics, and fusion models (with clinical data and radiomics features), using support vector machine (SVM) algorithm. The prediction performance of the models was assessed by accuracy, sensitivity, specificity, F1 score, and receiver operating characteristic (ROC) curve. Shapley Additive exPlanations (SHAP) was applied to analyze the interpretability and predictor importance of the model. RESULTS: A total of 1454 texture features were extracted from the NCCT images. In the test cohort, the ROC analysis showed that the radiomics model and the fusion model showed AUCs of 0.705 and 0.857, which outperformed the clinical model (0.643), with the fusion model exhibiting the best performance. Additionally, the accuracy and sensitivity of the fusion model were also the best among the models (84.8% and 93.8%, respectively). CONCLUSIONS: The model based on NCCT radiomics and machine learning has high predictive efficiency for the prognosis of AIS patients receiving conventional treatment, which can be used to assist early personalized clinical therapy.

Convolutional neural network for detecting rib fractures on chest radiographs: a feasibility study.

BACKGROUND: Chest radiography is the standard investigation for identifying rib fractures. The application of artificial intelligence (AI) for detecting rib fractures on chest radiographs is limited by image quality control and multilesion screening. To our knowledge, few studies have developed and verified the performance of an AI model for detecting rib fractures by using multi-center radiographs. And existing studies using chest radiographs for multiple rib fracture detection have used more complex and slower detection algorithms, so we aimed to create a multiple rib fracture detection model by using a convolutional neural network (CNN), based on multi-center and quality-normalised chest radiographs. METHODS: A total of 1080 radiographs with rib fractures were obtained and randomly divided into the training set (918 radiographs, 85%) and the testing set (162 radiographs, 15%). An object detection CNN, You Only Look Once v3 (YOLOv3), was adopted to build the detection model. Receiver operating characteristic (ROC) and free-response ROC (FROC) were used to evaluate the model's performance. A joint testing group of 162 radiographs with rib fractures and 233 radiographs without rib fractures was used as the internal testing set. Furthermore, an additional 201 radiographs, 121 with rib fractures and 80 without rib fractures, were independently validated to compare the CNN model performance with the diagnostic efficiency of radiologists. RESULTS: The sensitivity of the model in the training and testing sets was 92.0% and 91.1%, respectively, and the precision was 68.0% and 81.6%, respectively. FROC in the testing set showed that the sensitivity for whole-lesion detection reached 91.3% when the false-positive of each case was 0.56. In the joint testing group, the case-level accuracy, sensitivity, specificity, and area under the curve were 85.1%, 93.2%, 79.4%, and 0.92, respectively. At the fracture level and the case level in the independent validation set, the accuracy and sensitivity of the CNN model were always higher or close to radiologists' readings. CONCLUSIONS: The CNN model, based on YOLOv3, was sensitive for detecting rib fractures on chest radiographs and showed great potential in the preliminary screening of rib fractures, which indicated that CNN can help reduce missed diagnoses and relieve radiologists' workload. In this study, we developed and verified the performance of a novel CNN model for rib fracture detection by using radiography.

Integrative analysis of genome-wide association study and chromosomal enhancer maps identified brain region related pathways associated with ADHD.

Attention deficit/hyperactivity disorder (ADHD) is among the most common childhood onset psychiatric behavioral disorders, and the pathogenesis of ADHD is still unclear. Utilizing the latest genome wide association studies (GWAS) data and enhancer map, we explored the brain region related biological pathways associated with ADHD. The GWAS summary data of ADHD was driven from a published study, involving 20,183 ADHD cases and 35,191 healthy controls. The brain-related enhancer map was collected from ENCODE and Roadmap Epigenomics (ENCODE + Roadmap) including 489,581 enhancers. Firstly, the chromosomal enhancer maps of four brain regions were aligned with the ADHD GWAS summary data in order to obtain enhancer SNPs. Then the significant enhancers SNPs were subjected to the gene set enrichment analysis (GSEA) for identifying ADHD associated gene sets. A total of 866 pathways and 4 brain tissues were analyzed in this study. We detected several candidate genes for ADHD, such as AHI1, ALG2 and DNM1. We also detected several candidate biological pathways associated with ADHD, such as Reactome SEMA4D in semaphorin signaling and Reactome NCAM1 interactions. Our findings may provide a novel insight into the complex genetic mechanism of ADHD.

eQTLs Weighted Genetic Correlation Analysis Detected Brain Region Differences in Genetic Correlations for Complex Psychiatric Disorders.

BACKGROUND: Psychiatric disorders are usually caused by the dysfunction of various brain regions. Incorporating the genetic information of brain regions into correlation analysis can provide novel clues for pathogenetic and therapeutic studies of psychiatric disorders. METHODS: The latest genome-wide association study (GWAS) summary data of schizophrenia (SCZ), bipolar disorder (BIP), autism spectrum disorder (AUT), major depression disorder (MDD), and attention-deficit/hyperactivity disorder (ADHD) were obtained from the Psychiatric GWAS Consortium (PGC). The expression quantitative trait loci (eQTLs) datasets of 10 brain regions were driven from the genotype-tissue expression (GTEx) database. The PGC GWAS summaries were first weighted by the GTEx eQTLs summaries for each brain region. Linkage disequilibrium score regression was applied to the weighted GWAS summary data to detect genetic correlation for each pair of 5 disorders. RESULTS: Without considering brain region difference, significant genetic correlations were observed for BIP vs SCZ (P = 1.68 × 10-63), MDD vs SCZ (P = 5.08 × 10-45), ADHD vs MDD (P = 1.93 × 10-44), BIP vs MDD (P = 6.39 × 10-9), AUT vs SCZ (P = .0002), and ADHD vs SCZ (P = .0002). Utilizing brain region related eQTLs weighted LD score regression, different strengths of genetic correlations were observed within various brain regions for BIP vs SCZ, MDD vs SCZ, ADHD vs MDD, and SCZ vs ADHD. For example, the most significant genetic correlations were observed at anterior cingulate cortex (P = 1.13 × 10-34) for BIP vs SCZ. CONCLUSIONS: This study provides new clues for elucidating the mechanism of genetic correlations among various psychiatric disorders.

Integrating genome-wide association study, chromosomal enhancer maps and element-gene interaction networks detected brain regions related associations between elements and ADHD/IQ.

OBJECTIVE: To explore the associations between chemical elements and attention deficit hyperactivity disorder (ADHD)/intelligence quotient (IQ). METHODS: We applied elements related gene set enrichment analysis (ERGSEA) to explore the relationships between elements and ADHD/IQ. The GWAS dataset of ADHD was derived from the Psychiatric Genomics Consortium, involving 55,374 individuals. The GWAS dataset of IQ was derived from UK Biobank web-based measure (n = 17,862), UK Biobank touchscreen measure (n = 36,257), CHIC consortium (n = 12,441) and five additional cohorts (n = 11,748). Enhancer-gene datasets of eight brain tissues consist of 935 individuals. Utilizing the published GWAS summary and eight brain region-related chromosomal enhancer maps to obtain the SNP association testing signals. The element-gene interaction datasets of 21 elements were downloaded from the comparative toxicogenomics database (CTD). RESULTS: ERGSEA observed significant associations between 4 elements and ADHD, such as Al at Hippocampus Middle (P value = 0.040), As at Angular Gyrus (P value = 0.007) and Na at Hippocampus Middle (P value = 0.026). Additionally, ERGSEA identified that 5 elements were associated with IQ, mainly including Al at Dorsolateral Prefrontal Cortex (P value = 0.017), As at Dorsolateral Prefrontal Cortex (P value = 0.004) and Pb at Germinal Matrix (P value = 0.045). CONCLUSION: Our study results provide novel clues for understanding the associations between elements and ADHD/IQ. This study also illustrated the good performance of ERGSEA approach for complex diseases.

A Genomewide Integrative Analysis of GWAS and eQTLs Data Identifies Multiple Genes and Gene Sets Associated with Obesity.

To identify novel susceptibility genes and gene sets for obesity, we conducted a genomewide expression association analysis of obesity via integrating genomewide association study (GWAS) and expression quantitative trait loci (eQTLs) data. GWAS summary data of body mass index (BMI) and waist-to-hip ratio (WHR) was driven from a published study, totally involving 339,224 individuals. The eQTLs dataset (containing 927,753 eQTLs) was obtained from eQTLs meta-analysis of 5,311 subjects. Integrative analysis of GWAS and eQTLs data was conducted by SMR software. The SMR single gene analysis results were further subjected to gene set enrichment analysis (GSEA) for identifying obesity associated gene sets. A total of 13,311 annotated gene sets were analyzed in this study. SMR single gene analysis identified 20 BMI associated genes (TUFM, SPI1, APOB48R, etc.). Also 3 WHR associated genes were detected (CPEB4, WARS2, and L3MBTL3). The significant association between Chr16p11 and BMI was observed by GSEA (FDR adjusted p value = 0.040). The TGCTGCT, MIR-15A, MIR-16, MIR-15B, MIR-195, MIR-424, and MIR-497 (FDR adjusted p value = 0.049) gene set appeared to be linked with WHR. Our results provide novel clues for the genetic mechanism studies of obesity. This study also illustrated the good performance of SMR for susceptibility gene mapping.

Assessing the Associations of Blood Metabolites With Osteoporosis: A Mendelian Randomization Study.

Context: Osteoporosis is a metabolic bone disease. The effect of blood metabolites on the development of osteoporosis remains elusive. Objective: To explore the relationship between blood metabolites and osteoporosis. Design and Methods: We used 2286 unrelated white subjects for the discovery samples and 3143 unrelated white subjects from the Framingham Heart Study (FHS) for the replication samples. The bone mineral density (BMD) was measured using dual-energy X-ray absorptiometry. Genome-wide single nucleotide polymorphism (SNP) genotyping was performed using Affymetrix Human SNP Array 6.0 (for discovery samples) and Affymetrix SNP 500K and 50K array (for FHS replication samples). The SNP sets significantly associated with blood metabolites were obtained from a reported whole-genome sequencing study. For each subject, the genetic risk score of the metabolite was calculated from the genotype data of the metabolite-associated SNP sets. Pearson correlation analysis was conducted to evaluate the potential effect of blood metabolites on the variations in bone phenotypes; 10,000 permutations were conducted to calculate the empirical P value and false discovery rate. Results: We analyzed 481 blood metabolites. We identified multiple blood metabolites associated with hip BMD, such as 1,5-anhydroglucitol (Pdiscovery < 0.0001; Preplication = 0.0361), inosine (Pdiscovery = 0.0018; Preplication = 0.0256), theophylline (Pdiscovery = 0.0048; Preplication = 0.0433, gamma-glutamyl methionine (Pdiscovery = 0.0047; Preplication = 0.0471), 1-linoleoyl-2-arachidonoyl-GPC (18:2/20:4n6; Pdiscovery = 0.0018; Preplication = 0.0390), and X-12127 (Pdiscovery = 0.0002; Preplication = 0.0249). Conclusions: Our results suggest a modest effect of blood metabolites on the variations of BMD and identified several candidate blood metabolites for osteoporosis.

A Genome-wide Expression Association Analysis Identifies Genes and Pathways Associated with Amyotrophic Lateral Sclerosis.

Amyotrophic lateral sclerosis (ALS) is a neurodegenerative disease with strong genetic components. To identity novel risk variants for ALS, utilizing the latest genome-wide association studies (GWAS) and eQTL study data, we conducted a genome-wide expression association analysis by summary data-based Mendelian randomization (SMR) method. Summary data were derived from a large-scale GWAS of ALS, involving 12577 cases and 23475 controls. The eQTL annotation dataset included 923,021 cis-eQTL for 14,329 genes and 4732 trans-eQTL for 2612 genes. Genome-wide single gene expression association analysis was conducted by SMR software. To identify ALS-associated biological pathways, the SMR analysis results were further subjected to gene set enrichment analysis (GSEA). SMR single gene analysis identified one significant and four suggestive genes associated with ALS, including C9ORF72 (P value = 7.08 × 10-6), NT5C3L (P value = 1.33 × 10-5), GGNBP2 (P value = 1.81 × 10-5), ZNHIT3(P value = 2.94 × 10-5), and KIAA1600(P value = 9.97 × 10-5). GSEA identified 7 significant biological pathways, such as PEROXISOME (empirical P value = 0.006), GLYCOLYSIS_GLUCONEOGENESIS (empirical P value = 0.043), and ARACHIDONIC_ACID_ METABOLISM (empirical P value = 0.040). Our study provides novel clues for the genetic mechanism studies of ALS.

Integrating genome-wide association study and expression quantitative trait locus study identifies multiple genes and gene sets associated with schizophrenia.

Schizophrenia is a serious mental disease with high heritability. To better understand the genetic basis of schizophrenia, we conducted a large scale integrative analysis of genome-wide association study (GWAS) and expression quantitative trait loci (eQTLs) data. GWAS summary data was derived from a published GWAS of schizophrenia, containing 9394 schizophrenia patients and 12,462 healthy controls. The eQTLs dataset was obtained from an eQTLs meta-analysis of 5311 subjects, containing 923,021 cis-eQTLs for 14,329 genes and 4732 trans-eQTLs for 2612 genes. Genome-wide single gene expression association analysis was conducted by SMR software. The SMR analysis results were further subjected to gene set enrichment analysis (GSEA) to identify schizophrenia associated gene sets. SMR detected 49 genes significantly associated with schizophrenia. The top five significant genes were CRELD2 (p value=1.65×10-11), DIP2B (p value=3.94×10-11), ZDHHC18 (p value=1.52×10-10) and ZDHHC5 (p value=7.45×10-10), C11ORF75 (p value=3.70×10-9). GSEA identified 80 gene sets with p values <0.01. The top five significant gene sets were COWLING_MYCN_TARGETS (p value <0.001) and CHR16P11 (p value <0.001), ACTACCT_MIR196A_MIR196B (p value=0.002), CELLULAR_COMPONENT_DISASSEMBLY (p value=0.002) and GRAESSMANN_RESPONSE_TO_MC_AND_DOXORUBICIN_DN (p value=0.002). Our results provide useful information for revealing the genetic basis of schizophrenia.

GWAS summary-based pathway analysis correcting for the genetic confounding impact of environmental exposures.

Genome-wide association study (GWAS)-based pathway association analysis is a powerful approach for the genetic studies of human complex diseases. However, the genetic confounding effects of environment exposure-related genes can decrease the accuracy of GWAS-based pathway association analysis of target diseases. In this study, we developed a pathway association analysis approach, named Mendelian randomization-based pathway enrichment analysis (MRPEA), which was capable of correcting the genetic confounding effects of environmental exposures, using the GWAS summary data of environmental exposures. After analyzing the real GWAS summary data of cardiovascular disease and cigarette smoking, we observed significantly improved performance of MRPEA compared with traditional pathway association analysis (TPAA) without adjusting for environmental exposures. Further, simulation studies found that MRPEA generally outperformed TPAA under various scenarios. We hope that MRPEA could help to fill the gap of TPAA and identify novel causal pathways for complex diseases.

Integrating Genome-Wide Association and eQTLs Studies Identifies the Genes and Gene Sets Associated with Diabetes.

AIM: To identify novel candidate genes and gene sets for diabetes. METHODS: We performed an integrative analysis of genome-wide association studies (GWAS) and expression quantitative trait loci (eQTLs) data for diabetes. Summary data was driven from a large-scale GWAS of diabetes, totally involving 58,070 individuals. eQTLs dataset included 923,021 cis-eQTL for 14,329 genes and 4,732 trans-eQTL for 2,612 genes. Integrative analysis of GWAS and eQTLs data was conducted by summary data-based Mendelian randomization (SMR). To identify the gene sets associated with diabetes, the SMR single gene analysis results were further subjected to gene set enrichment analysis (GSEA). A total of 13,311 annotated gene sets were analyzed in this study. RESULTS: SMR analysis identified 6 genes significantly associated with fasting glucose, such as C11ORF10 (p value = 6.04 × 10-8), MRPL33 (p value = 1.24 × 10-7), and FADS1 (p value = 2.39 × 10-7). Gene set analysis identified HUANG_FOXA2_TARGETS_UP (false discovery rate = 0.047) associated with fasting glucose. CONCLUSION: Our study provides novel clues for clarifying the genetic mechanism of diabetes. This study also illustrated the good performance of SMR approach and extended it to gene set association analysis for complex diseases.

Genome-wide DNA methylation profile analysis identifies differentially methylated loci associated with ankylosis spondylitis.

BACKGROUND: Ankylosing spondylitis (AS) is a chronic rheumatic and autoimmune disease. Little is known about the potential role of DNA methylation in the pathogenesis of AS. This study was undertaken to explore the potential role of DNA methylation in the genetic mechanism of AS. METHODS: In this study, we compared the genome-wide DNA methylation profiles of peripheral blood mononuclear cells (PBMCs) between five AS patients and five healthy subjects, using the Illumina Infinium HumanMethylation450 BeadChip. Quantitative real-time reverse transcription-polymerase chain reaction (qRT-PCR) was performed to validate the relevance of the identified differentially methylated genes for AS, using another independent sample of five AS patients and five healthy subjects. RESULTS: Compared with healthy controls, we detected 1915 differentially methylated CpG sites mapped to 1214 genes. The HLA-DQB1 gene achieved the most significant signal (cg14323910, adjusted P = 1.84 × 10-6, ß difference = 0.5634) for AS. Additionally, the CpG site cg04777551 of HLA-DQB1 presented a suggestive association with AS (adjusted P = 1.46 × 10-3, ß difference = 0.3594). qRT-PCR observed that the mRNA expression level of HLA-DQB1 in AS PBMCs was significantly lower than that in healthy control PBMCs (ratio = 0.48 ± 0.10, P < 0.001). Gene Ontology (GO) and KEGG pathway enrichment analysis of differentially methylated genes identified four GO terms and 10 pathways for AS, functionally related to antigen dynamics and function. CONCLUSIONS: Our results demonstrated the altered DNA methylation profile of AS and implicated HLA-DQB1 in the development of AS.

Integrating genome-wide association study and expression quantitative trait loci data identifies multiple genes and gene set associated with neuroticism.

Neuroticism is a fundamental personality trait with significant genetic determinant. To identify novel susceptibility genes for neuroticism, we conducted an integrative analysis of genomic and transcriptomic data of genome wide association study (GWAS) and expression quantitative trait locus (eQTL) study. GWAS summary data was driven from published studies of neuroticism, totally involving 170,906 subjects. eQTL dataset containing 927,753 eQTLs were obtained from an eQTL meta-analysis of 5311 samples. Integrative analysis of GWAS and eQTL data was conducted by summary data-based Mendelian randomization (SMR) analysis software. To identify neuroticism associated gene sets, the SMR analysis results were further subjected to gene set enrichment analysis (GSEA). The gene set annotation dataset (containing 13,311 annotated gene sets) of GSEA Molecular Signatures Database was used. SMR single gene analysis identified 6 significant genes for neuroticism, including MSRA (p value=2.27×10-10), MGC57346 (p value=6.92×10-7), BLK (p value=1.01×10-6), XKR6 (p value=1.11×10-6), C17ORF69 (p value=1.12×10-6) and KIAA1267 (p value=4.00×10-6). Gene set enrichment analysis observed significant association for Chr8p23 gene set (false discovery rate=0.033). Our results provide novel clues for the genetic mechanism studies of neuroticism.

Tissue-specific pathway association analysis using genome-wide association study summaries.

MOTIVATION: Pathway association analysis has made great achievements in elucidating the genetic basis of human complex diseases. However, current pathway association analysis approaches fail to consider tissue-specificity. RESULTS: We developed a tissue-specific pathway interaction enrichment analysis algorithm (TPIEA). TPIEA was applied to two large Caucasian and Chinese genome-wide association study summary datasets of bone mineral density (BMD). TPIEA identified several significant pathways for BMD [false discovery rate (FDR) < 0.05], such as KEGG FOCAL ADHESION and KEGG AXON GUIDANCE, which had been demonstrated to be involved in the development of osteoporosis. We also compared the performance of TPIEA and classical pathway enrichment analysis, and TPIEA presented improved performance in recognizing disease relevant pathways. TPIEA may help to fill the gap of classic pathway association analysis approaches by considering tissue specificity. AVAILABILITY AND IMPLEMENTATION: The online web tool of TPIEA is available at https://sourceforge.net/projects/tpieav1/files CONTACT: fzhxjtu@mail.xjtu.edu.cnSupplementary information: Supplementary data are available at Bioinformatics online.

Comparative analysis of gene expression profiles of hip articular cartilage between non-traumatic necrosis and osteoarthritis.

Hip cartilage destruction is consistently observed in the non-traumatic osteonecrosis of femoral head (NOFH) and accelerates its bone necrosis. The molecular mechanism underlying the cartilage damage of NOFH remains elusive. In this study, we conducted a systematically comparative study of gene expression profiles between NOFH and osteoarthritis (OA). Hip articular cartilage specimens were collected from 12 NOFH patients and 12 controls with traumatic femoral neck fracture for microarray (n=4) and quantitative real-time PCR validation experiments (n=8). Gene expression profiling of articular cartilage was performed using Agilent Human 4×44K Microarray chip. The accuracy of microarray experiment was further validated by qRT-PCR. Gene expression results of OA hip cartilage were derived from previously published study. Significance Analysis of Microarrays (SAM) software was applied for identifying differently expressed genes. Gene ontology (GO) and pathway enrichment analysis were conducted by Gene Set Enrichment Analysis software and DAVID tool, respectively. Totally, 27 differently expressed genes were identified for NOFH. Comparing the gene expression profiles of NOFH cartilage and OA cartilage detected 8 common differently expressed genes, including COL5A1, OGN, ANGPTL4, CRIP1, NFIL3, METRNL, ID2 and STEAP1. GO comparative analysis identified 10 common significant GO terms, mainly implicated in apoptosis and development process. Pathway comparative analysis observed that ECM-receptor interaction pathway and focal adhesion pathway were enriched in the differently expressed genes of both NOFH and hip OA. In conclusion, we identified a set of differently expressed genes, GO and pathways for NOFH articular destruction, some of which were also involved in the hip OA. Our study results may help to reveal the pathogenetic similarities and differences of cartilage damage of NOFH and hip OA.

Lack of Association between the TSPAN18 Gene and Schizophrenia Based on New Data from Han Chinese and a Meta-Analysis.

Tetraspanin-18 (TSPAN18) potentially plays a role in the calcium signaling that is associated with dopamine-induced cortical neuron apoptosis and is considered to be an important mechanism in the pathogenesis of schizophrenia (SCZ). Furthermore, a genome-wide association study (GWAS) identified TSPAN18 as a possible susceptibility gene for SCZ. To validate these findings and reveal the effects of different inheritance models, seven single nucleotide polymorphisms (SNPs) of the TSPAN18 gene were analyzed in 443 patients with SCZ and 628 controls of Han Chinese descent via the SNPscan method. Single SNP, genotype, and association analyses with different models (i.e., additive, dominant, and recessive models) were performed, and the published datasets (2062 cases and 2053 controls) were combined with our results to determine the inheritance effects of the SNPs on SCZ. We observed genotypes and allele distributions of TSPAN18 gene did not show any significant associations in the Han Chinese population based on our experimental and meta-analytical results. Our findings indicate that the TSPAN18 gene is unlikely to be a major susceptibility gene for schizophrenia in Han Chinese.